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ABSTRACT 

The afterglow of a gamma-ray burst (GRB) is commonly thought to be due to con- 
tinuous deceleration of a relativistically expanding fireball in the surrounding medium. 
Assuming that the expansion of the fireball is adiabatic and that the density of the medium 
is a power-law function of shock radius, viz., next oc R~^, we analytically study the effects 
of the first-order radiative correction and the nonuniformity of the medium on a GRB 
afterglow. We first derive a new relation among the observed time, the shock radius and 
the fireball's Lorentz factor: ^0 = i?/4(4 — /c)7^c, and also derive a new relation among 
the comoving time, the shock radius and the fireball's Lorentz factor: t^o = 2/2/ (5 — k)'jc. 
We next study the evolution of the fireball by using the analytic solution of Blandford 
and McKee (1976). The radiation losses may not significantly influence this evolution. 
We further derive new scaling laws both between the X-ray flux and observed time and 
between the optical flux and observed time. We use these scaling laws to discuss the 
afterglows of GRB 970228 and GRB 970616, and find that if the spectral index of the 
electron distribution is p = 2.5, implied from the spectra of GRBs, the X-ray afterglow of 
GRB970616 is well fitted by assuming k = 2. 

Subject headings: gamma-ray : bursts — radiation mechanisms : nonthermal 
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1. INTRODUCTION 



The popular theoretical explanation for cosmological gamma-ray bursts (GRBs) is 
based on the fireball model. In this model, a GRB is thought to result from the dissipation 
of the kinetic energy of a relativistically expanding fireball. This dissipation can be either 
due to internal shocks formed during the collision between the shells with different Lorentz 
factors in the fireball (Rees & Mcszaros 1994; Paczynski & Xu 1994; Sari & Piran 1997), or 
due to external shocks formed by the fireball colliding with the surrounding medium (Rees 
& Meszaros 1992; Meszaros & Rees 1993; Katz 1994; Sari, Narayan, & Piran 1996). After 
the main GRB event, the expanding fireball is predicted to produce delayed emission at 
longer wavelengths (Paczynski & Rhoads 1993; Katz 1994; Meszaros & Rees 1997; Wijers, 
Rees, & Meszaros 1997; Reichart 1997; Waxman 1997a,b; Vietri 1997a,b; Tavani 1997; 
Sari 1997). It is easily understood that such an X, optical and/or radio afterglow is in 
fact due to continuous deceleration of the expanding fireball. 

Fortunately, the afterglows of GRBs have been detected recently in the error boxes 
at the sites of seven GRBs, e.g. GRB970228, GRB970402, GRB970508, GRB970616, 
GRB970815, GRB970828 and GRB971214. The current fireball model to explain these 
afterglows has in fact been divided into two sub-models. In the first sub-model, the 
fireball expansion following a GRB is thought to be adiabatic (Wijers et al. 1997; Reichart 
1997; Waxman 1997a,b). This is a reasonable assumption if the timescale for cooling of 
the accelerated electrons behind the shocks rapidly becomes longer than the expansion 
timescale of the fireball. The adiabatic expansion model has given a scahng relation 
between the expansion Lorentz factor and observed time: 7 oc t®^^^, provided that the 
surrounding medium is uniform and the expansion is ultra-relativistic. The studies of 
Wijers et al. (1997), Reichart (1997) and Waxman (1997a,b) showed that the adiabatic 
expansion model may satisfactorily explain the long-term behavior of the afterglows of 
GRB970228 and GRB970508. The more detailed numerical calculations by Huang et al. 
(1997) are consistent with these studies. An alternative sub- model accounted possibly 
for several properties of the afterglows of these two bursts has been presented by Vietri 
(1997b), who postulated that the expansion is highly radiative. This requires that the 
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accelerated electrons always cool more rapidly than the fireball expands. 

The purpose of the present work is to study analj^ically the effects of radiative cor- 
rections and nonuniformity of the surrounding matter on a GRB afterglow based on the 
adiabatic expansion model. Our study is stimulated by two motivations. First, Sari (1997) 
recently found that radiation losses may significantly influence the hydrodynamical evolu- 
tion of the fireball in the uniform medium. Second, it is possible that the sources of GRBs 
are merging of neutron star binaries (Narayan et al. 1992; Vietri 1996), failed supernovae 
(Woosley 1993), accretion-induced phase transitions of neutron stars (Cheng & Dai 1996), 
or hypernovae (Paczyhski 1998). If so, one may not think that the surrounding matter is 
uniform. A massive star as a progenitor of a GRB has possibly produced a stellar wind or 
a supernova remnant exists around a neutron star which is a source of the GRB. The stel- 
lar wind and/or the supernova remnant constitute the nonuniform surrounding medium 
of the GRB. It can be expected that the nonuniformity of the medium can shorten or 
prolong the relativistic expansion of the fireball. Such a behaviour should be compared 
with observations of afterglows of GRBs. 

In order to reach the purpose of the present work, in the next section we first derive a 
new relation among the observed time, the shock radius and the fireball's Lorentz factor, 
assuming that the density of the surrounding medium is a power-law function of shock 
radius. We also obtain a new relation among the comoving time, the shock radius and 
the fireball's Lorentz factor. We show that, in the case of uniform medium, the former 
relation turns out ot be consistent with the expression derived by Sari (1997), but the 
latter relation is different from the usual expression by a factor of 5/2. We next study 
the evolution of the relativistic fireball in the nonuniform medium by using the analytic 
solution found by Blandford & McKee (1976), and investigate the effect of the first-order 
radiative correction on the evolution. We find this effect may be insignificant. In section 
3, we derive the X-ray fiux from the fireball as a function of observed time, and give 
a scaling relation between the optical fiux and observed time. In section 4, we use our 
model to discuss the afterglows of GRB970228 and GRB970616. In the final section, we 
make a summary for our work. 
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2. THREE MEASURES OF TIME AND FIRST-ORDER RADIATIVE CORRECTION 

Even though the source of GRBs is unknown, one believes the following scenario for 
GRBs: A compact source (~ 10^ cm) releases an energy E comparable to that observed in 
gamma rays, E ~ 10^^ erg, over a time less than 100s. The huge optical depth in the source 
results in an initial fireball that expands and accelerates to relativistic velocity. During 
the acceleration, the fireball energy is converted to bulk kinetic energy with a Lorentz 
factor f] ~ 300. Of course, internal shocks may be formed in this period. Subsequently, 
the expansion of the fireball starts to be significantly influenced by the swept-up medium 
and two external shocks are formed: a forward blast wave and a reverse shock (Rees & 
Meszaros 1992). A GRB may be produced by nonthermal processes in shocks. After the 
GRB, the relativistic blast wave continues to sweep up the medium. The fireball may go 
into an adiabatic expansion phase. For simplicity, we assume that the medium density 
varies with shock radius based on the following expression: 



where Uq is the medium density at the radius Rq at which the fireball starts to be signif- 
icantly influenced by the medium. As usual, Ro is defined as 



where £^51 = E'/lO^^erg, 7^300 = r^/SOO, and no is in units of 1cm ^. 

The use of a power-law dependence of the density surrounding a GRB source is mo- 
tivated by the following points. First, some theoretical studies of supervova remnants 
(for a review see Ostriker & McKee 1988) suggest a power-law dependence of the rem- 
nant density at radius larger than some value. Second, the actual dependence might be 
more complicated than a power-law form. However, only in the case with a power-law 
form can we analytically study the evolution of a relativistic fireball. A study of the 
relatvistic-fireball expansion in the nonuniform medium without a power-law dependence 
of the matter density on radius needs numerical simulations and thus is beyond the scope 
of this paper. 




(1) 




(2) 
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2.1. Three Measures of Time 



We assume that the expansion is ultra-relativistic, and radiation losses are small. This 
implies that when most of the energy has been given to the medium, the energy in the 
shocked medium is constant and approximately equal to E. The rest mass of the shocked 
medium M oc R^"''. Since the medium was thermalized by the relativistic blast wave, 
its energy in the observer's frame is ~ M^'^ ~ E, where 7 is the Lorentz factor of the 
shocked medium just behind the shock. Thus, we have the following scaling law 

7 (X . (3) 

It should be emphasized that for a relativistic strong blast wave the shock's Lorentz factor 
r = \/27 (Blandford & McKee 1976). 

According to the scaling law (3), we next investigate relations among three different 
measures of time denoted as t, tco and which are measured in the rest frame of the 
burster, in the frame comoving with the fireball, and in the observer's frame respectively. 
When the shock propagates a small distance 5R ^ c5t, photons that are emitted from the 
shock will be observed on the observed timescale of 5t(^ = 5R/{2T'^c) = dR/{4'-f'^c) while 
the change in time in the frame comoving with the fireball is 6tco = 6R/{'yc). Integrating 
these equations over time and using the scaling law (3), we obtain 

R t , , 

^® " 4(4 - k)^^c " 4(4 - A;)72 ' ^ ' 

and 

- _ 2t 

~ (5 - k)^c ~ (5 - k)^ • ^ ^ 

In the case of uniform medium, viz., k = 0, eq. (4) turns out to be consistent with the 

expression derived by Sari (1997). In this case, comparing cq. (5) with the commonly 
used expression R/{'yc) (Waxman 1997a,b; Vietri 1997b), we see the present expression 
is a factor of 5/2 smaller. 

2.2. First-Order Radiative Correction 
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After having eqs. (4) and (5), we can study the adiabatic evolution of the fireball in 
the nonuniform medium, and investigate the effect of the first-order radiative correction 
on the evolution. 

The starting point of our study is eq. (69) of Blandford & McKee (1976) about the 
total energy of the fireball: 

167^?^ext"^pc27^-R^ 

^ = u^k ' 

which is constant (to lowest order). Here next is to be interpreted as the external density 
at the position of the shock. Combining eq. (6) with eqs. (1) and (4), we get 

7 = 7oh^ , (7) 



where 70 and to are defined as 



7o 

and 



{17 -Ak)E ^"^'-''^ 



(8) 



Accordingly, the shock radius can be written as 

/-^ X l/(4-fc) 

R = R,,l (I) . (10) 

It is easily seen that in the case of A; = eqs. (7) and (10) are in agreement with those 
derived by Sari (1997). 

In the following we consider only synchrotron emission from the accelerated electrons. 
We neglect the contribution of inverse- Compton (IC) emission from these electrons. This 
is becasue IC emission is not of importance particularly at late times of the fireball 
expansion (for a discussion to be seen in the final paragraph of this subsection) . In order 
to calculate the effect of synchrotron emission, we need to determine the magnetic field 
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strength and electron energy. As usual, we assume that the magnetic energy density in the 
comoving frame is a fraction of the total thermal energy density e' — A'j^ne^m.pC^ , viz., 
B' — (Stt^bc')^/^, and that the electrons carry a fraction of the energy. This implies 
that the Lorentz factor of the random motion of a typical electron in the comoving frame 
is 7em = iel^pl'tT^e- The ratio of the comoving-frame expansion time, t^o — 2i?/ (5 — A;)7C, 
to the synchrotron cooling time, tgyn = Girmec/ aT'yemB'^ , is 

2/(4-fe) 



tr 



t: 



syn 



where ti has been defined based on 



,2/(4-fc) 
''1 



128(4- A:) /m 



Ce^B{crTnocto)% 



,3-2ifc+2/(4-fe) 



3(5 - k) 

Thus, we can define the radiative efficiency as 



(12) 



syn 



(^lAe) 



2/(4-fe) 



''CO 



"syn 



(13) 



(tiAe)^/(^-'=) + 1 ■ 

This parameter accounts for a fraction of the energy of the accelerated electrons which 
is radiated away by the synchrotron emission. For <^ ti, ^ ^ 1; but for » ti, 
^ ^ {ti/t(sf'^'^-^\ Following eq. (84) of Blandford & McKee (1976), we derive the energy 
loss rate during the deceleration which is given by 4^(4 — k)'^7r^^cH'^nextmpC^ multiplied 
by CeC; where eq. (4) has been used. Using eqs. (7)-(9), we get the total power radiated 
per unit time: 

Ak E 



dE__ 21 

dt^' ^'^A{A-k)t^- 

Integrating this equation over time, we further obtain 



(14) 



(15) 



where tin is the observed initial time of the afterglow which is of the order of magnitude 
~ 10 s. We now define the parameter / as 



/ 



2/(4-fc) 



+ 1 



(16) 
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This parameter is in fact the ratio of the fireball energy at ^ ti to the initial total 
energy, and thus accounts for the effect of the first-order radiative correction. To estimate 
/, we adopt the following values: £'51 = 4, 77300 = 1, no = 1 cm~^, — 0.1, and ~ 0.1- 
0.3. Table 1 gives the values of 70, to, ti and / for different k. It can be seen from this 
table that for ~ 0.1 the effect of the radiative corrections may be insignificant, and 
for ~ 0.3 this effect may not be of importance in the case of A; > 0. This conclusion 
disagrees with that of Sari (1997), who didn't consider the factor ^ in eq. (14). For this 
reason, we will not take into account this effect in the next text. 

In the remainder of this subsection, we want to discuss the validity of two of our 
assumptions. First, the fireball expansion has been assumed to be adiabatic. As described 
in the introduction, this assumption in fact requires that the timescalc for cooling of the 
accelerated electrons behind the shocks rapidly becomes longer than the fireball expansion 
timescale. Now we assume that the initial expansion of the fireball is radiative. In this 
case (e.g., see Vietri 1997b), the relations among the shock radius, comoving-frame time 
and observed time are it! = 4(7 — 2/c)7^ci© and R— {A — k)^ctco, and the Lorentz factor 
of the fireball decreases as 

where 

i - - ^ pl/3 -1/3 -8/3 (.r.. 

Thus, the ratio of the comoving-frame expansion timescale to the synchrotron cooling 
timescale is given by 

T^-ir)^'^ (19) 



where tir is defined through the following expression 

34(7 -2A;) /m^ 
3(4 - k) W 



^(5-.)/(r-2.) ^ 64(7-2^ / rn,^ ^ . (20) 
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The values of tir for different k are calculated in Table 2. According to this table and eq. 
(19), we see that igyn > ^co for > ti,-, and therefore conclude that even if the fireball 
starts off with radiative dynamics, the transition to adiabatic evolution comes soon after 
the GRB. 

Second, we have not considered the contribution of IC emission to cooling of the ac- 
celerated electrons. It is well known that whether IC emission is important depends on 
the ratio of the IC power to synchrotron power: 

y-^, (21) 

where e'g and e'^ arc the synchrotron-photon and magnetic-field energy densities. For 
^co < ^syn and when emission is dominated by the synchrotron process, as argued by 
Waxman (1997a), the energy density e'g is a fraction tco/4tsyn of the electron energy 
density, and thus the ratio y is given by 

where t2 has been defined by 

/ ^ \ (4-fc)/2 

Therefore, we can conclude from Table 1 that IC emission is not an important process for 
coohng of the accelerated electrons behind the shocks for tf^ > t2. 

3. THE X-RAY AND OPTICAL RADIATION 

We first study the X-ray fiux from a relativistic fireball expanding in the nonuniform 
medium. The total power radiated per unit time has been given by eq. (14). In order 
to calculate the X-ray fiux, we need to consider radiation mechanisms of the accelerated 
electrons behind the shock. As shown in the last section, the synchrotron emission is the 
main mechanism for cooling of these electrons. We assume that the electron distribution 
behind the shock is power-law: 

OC , for 7em < 7e < 7e,max • (24) 



The estimate of 7e,max can be obtained by equating the electron acceleration timescale, 
ta — le'iTT'ec/ gB' , to the synchrotron cooling timescale so that 



7e,: 



max 



/ fiyre \ 



It should be pointed out that even though the maximum Lorentz factor estimated by this 
equation is rather high, the effect of IC emission on 7e,max can be neglected, as seen in 
eq. (22). For the electron distribution with eq. (24), the fraction of synchrotron power 
radiated in the X-ray region is (Vietri 1997a,b) 

.(3-p)/2 _ .{3-p)/2 



€}, ^"^ — e 

Cmax 

where e„ and ei are the upper and lower limits of the BeppoSAX instruments, 2 and lOkeV 
respectively, and emax is given by 

emax - ^7e,max7 = I6O7 MeV . (27) 

rUeC 

Inserting eq. (27) into eq. (26), we find 

f,^[6xio-yjf-^y\ (28) 

where the number in the brackets is two orders smaller than that of Vietri (1997b). 
Therefore, the expected X-ray flux is 

^'=wM (T^ri^f' ^''^ 

where d is the source distance, the constant has been computed for k = 0, p = 2.5, 
= 0.1, and the values of the other parameters used in Table 1. For <C in the 
above equation, 

3-k 3-p 

a — 1 — - — — — - — ; (30) 
8-2A; 2 ' ^ ^ 



but for ^0 ^ ti, 

6 — A; 3 — k 3 — p 
" ~ A-k ~ 8 -2k 2 ■ 
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(31) 



We next discuss the optical flux from the accelerated electrons behind the shock. For 
the electron distribution with eq. (24), the observed frequency of synchrotron emission 
at peak flux is 

„ - A^^2 ^ _ 2.1/2 .-A: fmpY einprnpcY^^ /'^©r'^' 

= S.O><.W(^)"(^)''^(^)-'%. (3.) 

where eqs. (8) and (9) have been used. Eq. (32) shows that is weakly dependent of 
k. The typical spectrum of the synchrotron emission has the following form: F^, oc i/", 
where a = 1/3 for z/ < z/^, and a = —{p — l)/2 for u > z/^. Since the comoving electron 
density n'g oc 7ncxt oc 

^-(3+.)/(8-2fc)^ and the comoving width of the fireball AT?' oc i?/7 oc 
^{5-k)/{8-2k) ^ then the comoving intensity II oc n'^B'AR' oc ^~(^+2'=)/(8-2fe) (^Meszaros & 
Rees 1997; Wijers et al. 1997). Thus, the observed peak flux as a function of time is 
Fi,^ oc t^'j^Il^ oc i^'^'^^^"^*^^ Once Um has entered the optical region, the observed flux at 
any frequency must vary according to 

F,^F,^{u/u^r^t-\ (33) 



where for u < Vm, 
2-k 



A-k 

for V > Urn, 



(34) 



It can be seen from eqs. (33)- (35) that only for A; < 2 the observed optical flux first 
increases and then decreases, but for k = 2 the flux is first kept to be constant and 
subsequently declines. 

4. DISCUSSION 
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We first use our model to discuss the X-ray abd optical afterglow of GRB970228. 
According to the observational results summarized by Wijers et al. (1997), we find p 2.4 
and /c by solving eqs. (31) and (35), which is well consistent with that of Waxman 
(1997a). The result of p 2.4 is consistent with the mean spectral index (p = 2.5) of 
GRBs measured in Band et al. (1993). This shows that the spectral index of the electron 
distribution due to shock acceleration is likely to be similar for GRBs and their afterglows. 

Our model can also be applied to discussing the X-ray afterglow of GRB970616. This 
burst was detected by BATSE on June 16.757UT. About 20min after the initial trigger, 
a transient X-ray XTE source was found in the error box of this burst (Connaughton et 
al. 1997), and 4 hours following the burst, scanning observations with the Proportional 
Counter Array on the Rossi X-ray Timing Explorer (RXTE) revealed an X-ray afterglow 
in the band 2-lOkeV with a flux ^ 1.1 x 10"" ergs"^ cm'^ (Marshall et al. 1997). On 
June 20.35UT, ASCA detected an X-ray flux from the XTE/IPN error box of GRB970616 
with ~ 3.7 X 10~^^ erg s~^ cm^^ in the band 0.7-7keV (Murakami et al. 1997). Using these 
values of the X-ray flux in eq. (29), we get a ^ 1.86. After knowing the value of a and 
assuming that p is equal to the mean spectral index of GRBs measured in Band et al. 
(1993), that is, p ~ 2.5, we solve eq. (35) and find k = 2. This result implies nonuniformity 
of the surrounding medium. 

We have found that the afterglows of GRB970616 and GRB970228 are weU explained 
by assuming two cases with k — 2 and respectively. We easily understand these two 
cases. In the first case, a neutron star as the GRB source has lain in a supernova remnant 
and/or a stellar wind due to the star's low velocity, so the postburst fireball has expanded 
in such a nonuniform medium; but in the second case, we conjecture that since the velocity 
of a neutron star as the progenitor of the GRB was very high, the GRB source has left a 
supernova remnant and/or a stellar wind, and the fireball has met the uniform interstellar 
medium. 

5. SUMMARY 

A GRB has been commonly believed to result from the dissipation of the kinetic energy 
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of a relativistically expanding fireball, and its X, optical and/or radio afterglow is due to 
continuous deceleration of the fireball. In this paper, we have assumed that the expansion 
of the fireball is adiabatic and ultra-relativistic. If compact objects (neutron stars or 
black holes) are an origin of the GRB, the surrounding medium of the fireball may be 
nonuniform due to the existence of a stellar wind and/or a supernova remnant. For 
simplicity, we have assumed that the density of the medium is a power-law function of 
shock radius, viz., next oc R~''. In addition, radiation losses may significantly influence 
the hydrodynamical evolution of the fireball (Sari 1997). In view of these two important 
arguments, we have analytically studied the effects of the first-order radiative correction 
and the nonuniformity of the medium on the GRB afterglow in this paper. The results 
of our study are summarized as follows: 

First, we have derived a new relation among the observed time, the shock radius and 
the fireball's Lorentz factor. We have also obtained a new relation among the comoving 
time, the shock radius and the fireball's Lorentz factor. We have shown that, in the case 
of uniform medium, the former relation turns out ot be consistent with the expression 
derived by Sari (1997), but the latter relation is smaller than the usually used expression 
by a factor of 5/2. 

Second, we have used the analytic solution of Blandford & McKee (1976) to derive 
the fireball's Lorentz factor and the shock radius as functions of observed time, which 
show that the nonuniformity of the medium must shorten {k < 0) or prolong {k > 0) 
the relativistic expansion of the fireball. Using these functions, we have further derived 

the radiation energy loss rate, and found the first-order radiative correction may be in- 
significant. This conclusion disagrees with that of Sari (1997) who neglected the radiative 
efficiency defined in eq. (13). 

Third, we have derived new scaling laws both between the X-ray fiux and observed 
time and between the optical flux and observed time. We have found that only for A; < 2 
the observed optical flux flrst increases and then decreases, but for k = 2 the flux is first 
kept to be constant and subsequently declines. 

Finally, We have used our model to discuss the afterglows of GRB 970616 and GRB 
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970228. We have seen that the afterglow of GRB970616 is well fitted by assuming k = 2. 
This value implies of the nonuniformity of the medium. 

We would like to thank the referee, Dr. Ralph Wijers, for very valuable suggestions, 
and Drs. K. S. Cheng and D. M. Wei for many helpful discussions. This work was 
supported by the National Natural Science Foundation of China. 
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Table 1. The typical values of some parameters for different k 
in the case of adiabatic evolution. 



To ^o(s) ^l,a(s) ti,b(s) /a /b ^2,a.(s) ^2,b(g) 

k = 4.67 3.31 X 10^ 200 1800 0.71 0.20 12.5 1013 

k = l 7.46 4.41 X 10^ 30 140 0.81 0.42 3.75 90.9 

k^ 2 18.56 6.61 x 10^ 5.0 15.0 0.94 0.73 1.25 11.3 



Notes — 7o, to, ti, f and t2 are defined by eqs. (8), (9), (12), (16) and (23) 
in the text, respectively. The values are computed for E^i = 4, tj^qq = 1, 
Uq = 1 cm^^ and = 0.1. The subscripts "a" and "b" represent 
= 0.1 and 0.3 respectively. 
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Table 2. The typical values of some parameters for different k 
in the case of radiative evolution. 







tor(s) 


tlr,a(s) 


tlr,b(s) 


k = 





0.227 


47.3 


220 


k = 


1 


0.317 


53.2 


210 


k^ 


2 


0.529 


47.9 


144 



Notes — tor and tir are defined by eqs. (18) and (20) in the text, 
respectively. The values of the parameters of the fireball and medium 
are taken as in Table 1. The subscripts "a" and "b" represent 
= 0.1 and 0.3 respectively. 
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